Alignment of galaxy spins in the vicinity of voids 
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We provide limits on the alignment of galaxy orientations witii the direction to the void center 
for galaxies lying near the edges of voids. We locate spherical voids in volume limited samples 
of galaxies from the Sloan Digital Sky Survey using the HB inspired void finder and investigate 
the orientation of (color selected) spiral galaxies that are nearly edge-on or face-on. In contrast 
with previous literature, we find no statistical evidence for departure from random orientations. 
Expressed in terms of the parameter c, introduced by Lee & Pen to describe the strength of such 
an alignment, we find that c < 0.11(0.13) at 95% (99.7%) confidence limit within a context of a toy 
model that assumes a perfectly spherical voids with sharp boundaries. 

PACS numbers; 98.80.Jk, 98.80.Cq 
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I. INTRODUCTION 

Within the currently accepted paradigm for structure 
formation, galaxies form through a complex series of 
mergers and accretion events in a beaded filamentary net- 
work of dark matter halos and subhalos. The cooling and 
condensation of gas within these dark matter hosts leads 
to the visible signatures by which we recognize galaxies 
P, [3, H, 0, 01 ■ While this broad picture has received im- 
pressive support in recent years, many of the detailed 
processes remain poorly understood. In particular the 
manner in which feedback alters the rate at which gas 
cools and the coupling between the angular momentum 
of the gas and the dark matter - both of which effect the 
sizes and shapes of disks - are outstanding problems in 
galaxy formation. We still cannot form realistic looking 
disk galaxies from an ab initio simulation in a cosmolog- 
ical context. 

This lack is important not just for our understanding 
of galaxy formation, but also because it impacts upon 
one of the most important probes of structure forma- 
tion and cosmology: weak gravitational lensing (see Q 
and references therein). A key assumption of most weak 
lensing analyses is that galaxies are randomly oriented. 
If the shapes of the observable parts of galaxies depend 
on their large-scale environment then this assumption is 
violated with potentially serious consequences for infer- 
ences derived from weak lensing methods. The sensitivity 
of weak lensing to violations of the random-orientation 
hypothesis is stronger if one attempts to use lensing to- 
mography 7] , as essentially all proposed surveys wish to 
do. 

There is agreement within the theoretical community 
that the shapes of dark matter halos do exhibit align- 
ment correlations, driven by the environment in which 
they form. This alignment is thought to be the theo- 
retical analog of the alignment correlation seen for mas- 
sive elliptical galaxies by d, HI, [l0|. Ell, and the mod- 
els found in are in good agreement with these lat- 



ter measurements. These alignments seem to persist to 
very large scales (^ 50 Mpc) which may indicate that 
the commonly adopted approximation of Gaussian tidal 
fields with negligible non-linear effects might not be fully 
justified [11, m. 

Both the observational case for, and the theoretical 
models of, the alignment of disks or angular momentum 
directions is less settled. Standard assumptions are that 
galactic disks lie perpendicular to the angular momentum 
direction of the baryons, or the inner halo, or the entire 
halo. Some hydrodynamic simulations show that that 
there is a strong, but not perfect, correlation between 
the angular momentum vector of gas and that of the dark 
matter halos hosting the galaxy 17:]. Other simulations 
indicate that while the disk is very well aligned with the 
halo orientation in the inner regions of the dark mater ha- 
los (at distances less than 10% of the virial radius), the 
angular momenta between the inner and outer regions are 
essentially uncorrelated [l^. In any case, it is reasonable 
to assume that the degree of alignment of galaxies with 
large-scale structure is smaller than that between dark 
matter halo angular momenta and the structure. One 
way of searching for such a signature is to look at the ori- 
entation of angular momentum vectors on the outskirts of 
cosmic voids. The gravitational evolution of matter fluc- 
tuations in the cosmological context tends to make under- 
densities (voids) rounder, while evolving over-densities 
(halos) into more elongated structures [l^. This larger 
degree of symmetry in the voids makes more amenable 
to analysis. The theoretical situation here is also com- 
plex, though iV-body simulations, in general, find a very 
small or negligible correlation between the angular mo- 
mentum of the inner parts of dark matter halos in thick 
shells around voids and the direction to the center of the 
void. In particular, the authors of references [l3| and [l9| 
found no alignment. The authors of reference J20| found 
an alignment only in very thin shells around the void for 
halos with quiescent merging histories, not for the thicker 
shells used in most observational work. They also claim 
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very strong numerical requirements for the method to 
converge. On the other hand, the authors of [2lj found 
an ahgnment if they considered the outer parts of the halo 
in their calculation of the angular momentum (but not 
the inner regions alone), even with thick shells and less 
stringent numerical requirements than [20j. Reference 
[23| does not study the same problem, but those authors 
found that halos in sheets tend to have angular momenta 
parallel to the sheets, though the overall level of align- 
ment is low. In [2^ the authors find a range of behaviours 
ranging from angular momentum pointing in the direc- 
tion perpendicular to the structure, to no-alignment, and 
then to a regime where the angular momentum is con- 
tained by the plane defined by the surrounding structure, 
depending on the halo mass. 

Given this theoretical background, the observational 
situation is intriguiging. Recently [24| presented a mea- 
surement from the 2dFGRS and SDSS that spiral galax- 
ies located on the shells of the largest cosmic voids have 
rotation axes that lie preferentially in the void surface. 
The strength of this alignment was much higher than all 
measurements of angular momentum correlations of dark 
matter halos in A'^-body simulations, which as we argued 
above are expected to be an upper limit to the corre- 
lations of galaxies. In this paper we try to update this 
analysis with the most recent SDSS data. 



where L = L/|L| and T is defined by 



TkiTki 
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with Tij = Tij — Tr(T)(5y /3. The parameter c which 
can take values between and 1 describes the effect of 
non-linearities. The value of c = 3/5 corresponds to T 
and I being fully uncorrelated. This ansatz is derived 
from "marginalising" over all possible moments of iner- 
tia in the limit that the tidal field and moment of iner- 
tia are completely uncorrelated. The parameter c was 
then introduced to account for the randomizing effects of 
small-scale physics. 

We can now build a simple toy model for the behaviour 
of angular momenta in the vicinity of voids within the 
context of TTT and ansatz of Equation fl]). Assuming 
a sharp boundary, the potential just outside the void is 
(j) (X +l/r in comoving coordinates. Taking derivatives, 
evaluating the resulting at (0,0, 1) and rescahng ac- 
cording to the Equation [5l one obtains 



f = diag(-l/%/6, -l/%/6, +2/V6). 



(6) 



This result now holds for any T in the frame of its eigen- 
vectors, with z pointing in the radial direction. Using 
Equation ^ we then find 



II. THEORY 

The most commonly invoked theory for the origin of 
galaxy spins is known at the Tidal Torque Theory (TTT) 
[25l . Ill, [13, m, nil, which essentially states that the 
dark matter spins up as a result of coupling of the local 
quadrupole moment of matter distribution to the shear 
field, giving 



Lj(t) = a^{t)beijkTjiIi 



Ik, 



(1) 



where a is the scale factor of the Universe, D is the 
growth factor and tiju is the Levi-Civita symbol. The 
local inertia tensor /y of the protohalo (the mass that 
will later form the dark matter halo) in Lagrangian space 
is given by 



3o / qiQjA^q, 
Jv 



(2) 



where qi are the Lagrangian coordinates around the cen- 
tre of mass of the halo and po is the mean density. The 
local shear tensor Tij is defined by 



(3) 



where (j) is the gravitational potential. 

Lee and Pen fsdl have proposed the following simple 
ansatz for the relation between the galaxy spin vector L 
and the shear tensor T 
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Following the literature, we next assume that the prob- 
ability distribution function for L can be described as a 
Gaussian, with zero mean and second moments {LiLj): 



P{L) 



1 



(27r)^/^ detC 



■ exp 



LC-^L 



(8) 



Writing L in spherical coordinates and integrating over 
radial and azimuthal coordinates, one can obtain the 
probability distribution function for P{0), where is the 
remaining zenith coordinate, namely the angle between 
the axis of L and radial direction from the centre of void: 



P{0) = sine 



2(l-c)V2 + c 



(2-|-3ccos2 6'-2c 



,3/2- 



(9) 



In this equation, P{9) is normalised so that the proba- 
bility integrates to one over the interval 9 — 0. . .7r/2. 
For c = 0, P{9) = sinO, which is a purely geometrical 
factor. A result applicable to a more general settings can 
be found in [31]. 

The purpose of this paper is to test whether there is 
any evidence that P{9) deviates from the simple sine form 
in the vicinity of voids and if so, what values can the as- 
sociated parameter c take. We stress again that the voids 
boundaries are not sharp (as seen in the Figure [5]) , but 
fitting the Equation © to the data should nevertheless 
give some idea about the size of the parameter c. 
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III. DATA AND METHOD 

200 - 

Our method is essentially the same as that of [l^l and 
is composed of three steps. First we identify the voids in 
the SDSS data. Second, we locate spiral galaxies in the 
vicinity of these voids. Third, we investigate if the distri- 
bution of spin axes of nearly face-on or edge-on galaxies 
is inconsistent with random distribution. These steps j-j 
are discussed in more details in the following subsec- 2 
tions. Throughout we use concordant cosmology with g 100 - 

r^A — 0.75, rim — 0.25 and the reduced Hubblc's con- « 
stant h = iJo/lOOkm/s/Mpc = 0.7. Our data come from 
the 6th data release of the SDSS [33| and we use the 
Value Added Galaxy Catalog (VAGC) [3J| to obtain the 59 - 

volume-limited samples that are used to find voids. We 
limit ourselves to the northern galactic cap. 



A. Void Finder 

The void finder is conceptually based on the HB void 
finder of [l^. We start by creating a volume limited 
catalogue of galaxies. We create two volume limited cat- 
alogues, depending on the limiting absolute extinction- 
corrected r-band Petrosian magnitude, which we choose 
to be either rri]^ — —22 or 111^ ~ —21. We use only 
parts of the survey where the survey completeness de- 
scribed by the SDSS parameter FGOTMAIN is above 0.82 
and randomly remove galaxies to make the sample uni- 
formly complete at this level. In addition we also use the 
VAGC's random catalogues to create random catalogues 
with the same selection function as the main catalogues, 
but with about ten times as many points. We then find 
voids using the following algorithm. 

1. Pick a random point in the survey volume. Lo- 
cate four closest neighbours in the real catalogue. 
Determine the sphere that is defined by these four 
points. 

2. Find if any other real galaxies are inside the sphere. 
If so, diminish the radius accordingly. 

3. If the sphere's radius is less than 10 Mpc/h for the 
niL = —21 or 14 Mpc/h for the ttil = —22, discard 
it and go back to 1. 

4. Count the number of random points inside the 
sphere and if this number is more than two sigma 
below the number expected given the number den- 
sity of random points, discard it and go back to 1. 
Moreover, if the first moment of radial vectors of 
random points inside the void is inconsistent with 
random, discard them and go back to 1. This step 
essentially ensures that the void is fully within the 
survey. 

5. Compare the void candidate to the existing voids. 
If it touches any of the existing voids, then dis- 
card the smaller of the two. Check if the new void 
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FIG. 1: This figure shows the distribution of voids sizes for 
rriL = —21 (top) and tul ~ —22 (bottom). Note that the 
minimum void size is an input parameter for the void finder 
and equals to 10 and 14 Mpc/h respectively. 

touches any other voids and discard those as well. 
This ensures that two close small voids are merged 
into a bigger one if possible. 

This process is repeated for all random points. Its re- 
sult is a catalogue of non-overlapping spherical voids that 
contain no galaxies that are brighter than ttil (though 
less bright galaxies can be in the voids). Although not 
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FIG. 2: This figure shows how the mean number density of 
galaxies change as a function of distance from the void cen- 
tre scaled by the void radius for the m_t ~ —21 (top) and 
niL = —22 (bottom) samples. Error-bars are derived from 
Monte Carlo simulations in which void centres are randomly 
positioned within the volume limited catalogue and are ap- 
proximate and heavily correlated. 



strictly ensured by the algorithm, we find that the order 
in which random points are taken does not matter - the 
number of void solutions that are affected by this is less 
than 0.5%. Some basic properties about the voids the we 
find are tabulated in the Table HI 



To aid systematic checks we also form an ensemble 
of random void catalogues. In this catalogue we take 
the real number of voids and their radii and distribute 
them non-overlappingiy, but otherwise randomly within 
the survey volume. 

Throughout we ignore the distinction between real 
space and redshift space. Since our voids are relatively 
large and separated from high-density peaks, this should 
be a relatively safe approximation. 

We plot basic properties of our voids in Figures [1] and 
[5] These are here to illustrate the properties of our voids 
in broad brushes and to show the dependence of voids on 
the void-finding algorithm properties. The main thrust 
of this paper is to analyze the properties of galaxy spins 
in the vicinity of voids; we will return to properties of 
voids in a forthcoming publication. 



B. Galaxies in the vicinity of voids 

We next locate spiral galaxies that reside in the vicinity 
of voids. We now again work with the full SDSS cata- 
logue and select spiral galaxies by a colour cut g—r < 0.6, 
where g and r are model magnitudes from the SDSS 
pipeline. Spiral galaxy is considered to be at the edge 
of a void if its distance to the void centre is between rvoid 
and rvoid + ^aper, where rtaper = 4Mpc/h is the taper 
length. We limit ourselves to either face-on or edge-on 
galaxies. To determine the inclination of the galaxy, we 
use the adaptive moments e+ and ex from the SDSS 
photometric data reduction and calculate the axis ratio 
q using [35]: 



1 - e 

1 + e 



1/2 



(10) 



where e ^ yej, + e^. Galaxies whose q > 0.96 are 

deemed to be face on and those with q < 0.27 are edge 
one. These criteria are somewhat more relaxed than 
those of [iSli but visual inspection of a random sub- 
sample ensured that we are indeed selecting face-on and 
edge-on spiral galaxies. For face-on galaxies, we assumed 
that the angular momentum axis is aligned with the ra- 
dial direction toward the galaxy. Similarly, for the edge- 
on galaxies it was assumed to lie in the plane perpendicu- 
lar to the radial direction and along the galaxies apparent 
minor axis. The latter was again inferred from the SDSS 
photometric reduction using the isophotal angle in the 
r-band and again visually confirmed to be measuring the 
correct quantity. 

For the sample of galaxies and voids we then measure 
the probability distribution function for 9, the angle be- 
tween galaxy spin and the radial vector from the centre of 
the void toward the galaxy. Since the axis of the galaxy's 
spin vector is a spin-2 quantity, the angle is constrained 
to lie between 0° and 90°. We therefore measure the 
probability distribution function for this angle in 20 bins 
of uniform size in this range. The distribution of 6 should 
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m.L = -21 


JUL = -22 


Number of voids 




438 






447 




Minimum void radius 




10 Mpc//i 






14 Mpc/fe 




Maximum void radius 




17.8 Mpc//i 






21.4 Mpc/h 




redshift range 




0.169 - 0.183 






0.111 - 0.123 






face-on 


edge-on 


both 


face-on 


edge-on 


both 


Number of galaxies 


255 


323 


578 


151 


107 


258 


NaiVe 


25.0 


16.4 


17.6 


21.8 


20.0 


22.5 


MC derived 


23.9 


14.1 


16.8 


20.0 


18.5 


21.8 


Constraints on c 


c < 0.11(0.13) 


c < 0.16(0.19) 



TABLE I: This table shows the basic properties of voids that we find in our data and the resulting values for null model. 
The null model has no degrees of freedom and there are 20 datapoints. The 95% limits correspond to of 9.6 - 34.2. The 
bottom section shows 95% (99.7%) limits on the c parameter. See text for discussion. 




FIG. 3: This figure shows the measured probability distribution for the angle 6 in the tul ~ —21 sample. Panels show results 
corresponding to face-on only (top), edge-on only (middle) and combined (bottom) spiral galaxies. The dashed black line 
corresponds to the theoretical expectation p{9) = sin6, while the solid red-line corresponds to the Monte Carlo results. The 
jaggedness in the red line is real and not a numerical artifact. Points are measurements with Monte Carlo derived errorbars. 
Dotted green lines are theoretical predictions for the models with c = 0.1, 0.3, 0.5. 
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FIG. 4: Same as Figure [3] but for a catalogue of voids derived from mi = —22 volume limited sample. 



follow sin 9 is our sampling of galaxies around voids would 
be random and if there were no edge effects associated 
with voids. However, in practice, while the voids were re- 
quired to lie within the survey, the voids fattened by the 
taper radius might lie outside survey. Moreover, the cos- 
mic web makes the sampling of positions non-uniform. In 
order to correct for these effects, we Monte-Carlo (MC) 
our errorbars as follows. For all galaxies that lie in the 
vicinity of void, including those that are neither face-on 
or edge-on, we randomly permute their spin vector prop- 
erties and calculate the resulting p{9) for an ensemble 
of permutations. The question we are therefore asking 
is: given the number of face-on, edge-on and remaining 
galaxies in the vicinity of voids, what is the probability 
that this distribution is not random? Using Monte-Carlo 
simulations we calculate the number of galaxies expected 
in each 9 bin and the correlations between adjacent bins. 
Moreover, under assumption of Gaussianity, this allows 
us to calculate the exact by inverting the correlation 
matrix Hol. 



IV. RESULTS 

Our main results are plotted in Figures [3] and U) These 
figures show the expected signal in red solid line. This 
was derived using Monte Carlo simulations explained in 
the previous section. The jaggedness in the line is real 
and corresponds to the concrete distribution of cosmic 
web around the particular voids that we find. The plot- 
ted errorbars are also derived from the Monte Carlo Sim- 
ulation. 

It is clear already from the picture, that there is no 
strong evidence for departure from expected signal. We 
quantify this in Table HI In this table we quote two differ- 
ent values. The nai've are calculated assuming that 
the number of galaxies in each bin is Poisson distributed 
and that bins are independent. The Monte-carlo derived 
X^s use the mean values for null theory derived from the 
Monte Carlo simulations and evaluate taking in ac- 
count the full covariance properties from MC samples. 
Since there are twenty data-points at which we estimate 
p{6), we note that there are no detections at greater than 
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2 sigma. 

We note that the correctly calculated is in fact sur- 
prisingly close to the naively calculated one. In fact, this 
is due to the large number of voids that can be found us- 
ing the DR6 dataset. If we artificially limit the number 
of voids, the difference can be much larger. For example, 
taking the first 50 voids of the tol = —22 sample, the 
for face-on galaxies drops from 45.3 using naive to 
17.9 using MC x^ with twenty degrees of freedom. This 
corresponds to a drop from over 3 sigma detection to less 
than 1 sigma. It is therefore crucial for the spinologist's 
peace of mind that x^s are calculated robustly. 

Next we calculate limits on the parameter c, which are 
also shown in the Table [D Two magnitude bins give very 
strong, but also very similar constraints on c. 



A. Systematics 

It is possible that there is some crucial aspect of the 
analysis which is adversely affecting our results. We 
therefore compare the stability of our results with re- 
spect to the following tweaks to our data reduction and 
find no significant change in our results: 

1. Including galaxies inside the void. 

2. Considering just galaxies inside the void. 

3. Changing the major to minor axes ratios required 
for face-on and edge-on criteria. In p articular, we 
tested using the same q ratios as [2J|, namely q > 
0.978 for face-on and q < 0.208 for edge-on. 

4. Removing the color cut. This significantly increases 
our catalogue of galaxies, which becomes heavily 
contaminated by ellipticals. 



reality, but the numbers that we report should give an 
idea about the size of allowed values of c. Our limits are 
consistent with most A^-body simulations, but there is 
some tension with the largest signal found by 21| , which 
is c~ 0.15. 

Our results are also inconsistent with previous results 
of 0, who find c = O.lto l- The reasons for this dis- 
crepancy are not immediately clear, but we note that 
the catalogue that we use is considerably larger and has 
a much better filling factor that dramatically increases 
the number of voids. It possible that the results of [l^l 
were simply a statistical fluctuation. 

When compared to the iV-body simulations, the tidal 
torque theory seems to give qualitatively, but not quan- 
titatively correct results [3^, |37| . The ansatz of Equation 
dU is perhaps less universally applicable. It was found 
to work reasonably well in early data [ssj : however, there 
was recently a tentative detection of correlation of spin 
vectors that should vanish in this theory [3i|. In this 
paper we find c is consistent with zero. Moreover, there 
does not seem to be good evidence for decoupling of scales 
between the tidal field and the local moment of inertia; 
a simplistic fc-counting gives the same scale dependence. 
Finally, it is not clear whether the parameter c should 
have one universal value or be an environment or scale- 
dependent variable. 

On the basis of our revisiting the alignment of galaxy 
spins in the vicinity of voids we are again able to reconcile 
results from A^-body simulations with observations. This 
is especially important in the light of the forth-coming 
weak lensing surveys; a nagging discrepancy between sim- 
ulations and observations might lower our trust in cor- 
rections and systematic errors due to intrinsic alignment 
that are derived from A^-body simulations. While more 
work must be done, the outlook is currently quite posi- 
tive. 



V. DISCUSSION & CONCLUSION 



In this paper we have reconsidered the alignment of 
galaxy spin axes in the vicinity of voids. We find no 
statistical departure from random for the orientation of 
galaxy spin axis in the vicinity of voids. Expressed in 
terms of the c parameter, we find very stringent limit on 
c < 0.11 (0.13) at 95% (99.7%) confidence limits. These 
limits assume a toy model with perfectly spherical voids 
and sharp edges. These are only approximately true in 
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